Minimum spanning trees on weighted scale-free networks 
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A complete understanding of real networks requires us to understand the consequences of the 
uneven interaction strengths between a system's components. Here we use the minimum spanning 
tree (MST) to explore the effect of weight assignment and network topology on the organization of 
complex networks. We find that if the weight distribution is correlated with the network topology, 
the MSTs are either scale-free or exponential. In contrast, when the correlations between weights 
and topology are absent, the MST degree distribution is a power-law and independent of the weight 
distribution. These results offer a systematic way to explore the impact of weak links on the structure 
and integrity of complex networks. 

PACS numbers: 89.75.Hc,89.75.Fb,89.75.Da,05.10.-a 
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The study of many complex systems have benefited 
from representing them as networks 0] , examples includ- 
ing metabolic networks Q, describing the reactions in 
a cell's metabolism; the protein interaction network 
capturing the binding interactions between a cell's pro- 
teins; and the World Wide Web and email networks 0,0 
linking web-pages or people together via URLs or emails. 
For these systems there is extensive empirical evidence 
indicating that the degree (or connectivity) distribution 
of the nodes follows a power-law, strongly influencing ev- 
erything from network robustness JQ to disease spreading 
pj. However, to fully characterize these systems, we need 
to acknowledge the fact that the links can differ in their 
strength and importance [1, El El El ■ Indeed, in a 
social network the strength of the relationship between 
two long-time friends differs from that between two ca- 
sual business associates ^3; m ecological systems the 
strength of a particular pair-interaction between species 
is crucial for population dynamics [l4| , ecosystem stabil- 
ity [lij and development in stressed environments |l6j . 
Thus in most networks the links are not binary (present 
or absent), but have a strength that quantifies the im- 
portance of the particular node-to-node interaction. 

The weakest links can carry particular significance in 
some weighted networks 13J. For instance, the speed of 
data transmission between two computers is limited by 
the link with the smallest bandwidth ("bottleneck"), or 
the activity of a metabolic pathway is determined by the 
rate of the slowest reaction. Furthermore, weak links can 
affect the overall network integrity. For example, ecolog- 
ical communities may experience dramatic effects upon 
the removal of weak interactors [l^. To systematically 
uncover the location and the role of weak links in a com- 
plex network, we use the minimum spanning tree (MST) , 
which for an N node network represents the loopless sub- 
graph of (N — 1) links that reaches all nodes while min- 
imizing the sum of the link weights [l7i fisl flfl Etl ] . By 
avoiding the strong links and preferentially following the 
weakest ones, the MST selects the lowest weight back- 



bone of a network. 

We start by examining the correlations between 
weights and network structure for several real systems, 
allowing us to construct a model system whose weight 
distribution mimics the statistical features of real net- 
works. We then show that the large-scale structure of 
the MSTs depends on the way the weights are placed 
in the network: For systems whose weight distribution 
is correlated with the network topology, the MSTs are 
either scale-free or exponential. In contrast, when the 
correlations between weights and topology are removed, 
the MST degree distribution is a power-law with a degree 
exponent close to the degree exponent of the original net- 
work, independent of the weight distribution. 

Topology Correlated Weights. — To uncover the func- 
tional relationship between network topology and link 
weights, in Fig. we display the dependence of the 
weights on the node degrees for the E. coli metabolic 
network, where the link weights represent the optimal 
metabolic fluxes H3; the US Airport Network (US AN) 
where the weights reflect the total number of passen- 
gers travelling between two airports between 1992 and 
2002; and the link betweenness-centrality (BC), repre- 
senting the number of shortest paths along a link for the 
Barabasi- Albert (BA) scale-free model |22j. For each of 
these systems the weight distributions follow a power- 
law |3| (not shown) and, as Fig^shows, the average link 
weight scales with the degrees of the nodes on the two 
ends of a link as (wy) ~ (kikj) e , similar to the scaling 
found for the World Airport Network 

These empirical observations allow us to assign weights 
to the links of a network for which we have only the 
network topology. To systematically study the role of 
the weight distribution on the structure of the MST 
we use several weight assignments, (i) First, we choose 
Wij = kikj (see inset Fig. [I). Note that the MST gen- 
erated by this weight assignment is identical to the MST 
obtained for weights w'^ — (u)y) e with an y 6 > 0, as it is 
the rank of the weights and not their absolute value that 
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FIG. 1: The average weight of a link between nodes i and 
j shown as function of the link end-point degree product 
kikj. The symbols represent (i) the US AN with the num- 
ber of passengers as link weights (filled squares); (ii) the E. 
coli metabolic network with optimized flux as link weight 
p?H (filled triangles); (iii) Barabasi- Albert scale-free model 
with betweenness-centrality (BC) as link weights (open cir- 



cles). The solid line (w ~ (kikj 



and the dashed line 



(w ~ (kikj) ) serve as guides to the eye. Inset: The 
weights are determined by the end-point degrees ki and kj\ 
(i) Wij — kikj, (ii) Wij = max(ki,kj), (iii) Wij = min(ki,kj) 
or the inverse thereof (iv)-(vi). 



determines the MST [23j. We have also studied the two 
extreme cases of topology-correlated weights, distributed 



according to (ii) Wi 



and (iii) Wi 



l , where 



k m in = min(fcj,fcj) and k max — max(fc;, kj) and with 
fc^ m < kikj < k^ lax . Finally, we investigated the struc- 
ture of the maximal spanning trees 24] for the above 
weight choices by determining the MST after transform- 
ing cases (i)-(iii) as u^. = 1/wy, resulting in the link- 
weight choices (iv) ~ 1 /kikj, (v) Wij ~ l/k max and 
(vi) Wij ~ 1/k 

nun • 

Weight Distributions. — To characterize the obtained 
weighted networks we first study their weight distribu- 
tion. For this, we grow scale-free networks according to 
the BA model 22], the resulting networks having a de- 
gree distribution P(k) ~ k^ 1 with 7 = 3. We then assign 
a weight to each link according to (i)-(vi). For networks 
whose degrees at the two ends of a link are uncorrelated 
we can determine the weight distribution analytically us- 
ing order statistics |25j , finding 
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FIG. 2: Distribution of link weights on N = 10° node scale- 
free networks. Link- weight choice (i) (triangles), (ii) (squares) 
and (iii) (circles) are all heavy tailed. The analytical predic- 
tions (Eqs. m - £3) are indicated as solid lines. Note that the 
solid curves have been shifted vertically without changing the 
character of the scaling law. Inset: The inverse weight distri- 
butions (iv)-(vi) (triangles, squares and circles respectively) 
and the analytical predictions shown as continuous lines. 



lations are obtained after the variable change w' = l/w. 
For Wij ~ (kikj) 6 the relationship between the exponent 



of the weight distribution (P(w) ~ w ",w ^> 1) 
exponent of the degree distribution 7 and 6 is 
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Corresponding expressions for the inverse degree corre- 



valid for cases (i) and (ii). In Fig. [3 we compare the nu- 
merically determined weight distributions with the scal- 
ing predicted by our analytical expressions, finding that 
the numerical curves display a w-dependency close to 
that of Eqs. CQ~ ©7 unaffected by the degree-degree 
correlations in the model [2H. \F\ . Note that power-law 
weight distributions like those in Fig. |2 have been ob- 
served for a wide range of network based dynamical pro- 
cesses [28). 

Minimum Spanning Trees. — The MSTs were generated 
using Prim's greedy algorithm |29| : starting from a ran- 
domly selected node, at each time step we add the link 
(and hence a node) with the smallest weight among the 
links connected to the already accepted nodes. When- 
ever m links with the same (smallest) weight are encoun- 
tered, we break the degeneracy by randomly selecting one 
among them with probability 1/m. 

The numerical results indicate that the degree distri- 
bution of the resulting MSTs fall into two distinct classes 
[30I l3lT | . Weight choices (i) and (ii) give rise to exponen- 
tial MST degree distributions (Fig. |3Jl), while choices 
(iii)-(vi) result in power-law distributed MST degrees 
(Fig. 1^)). We can understand the exponential nature 
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FIG. 3: Degree and weight distribution of TV = 10 node 
MSTs. (a) The degree distribution for weights proportional 
to either kikj (i) or k ma x (ii) are dominated by an exponential 
cut-off, while (b) it is heavy-tailed for weights proportional 
to kmi n (iii) and inversely proportional to either kikj (iv), 
kmax (v) and k m i n (vi). (c) The distribution of link weights 
on the MSTs is a power law for (i), (ii) and (vi) (w x 10 3 ), 
while (d) it is dominated by an exponential cut-off for (iii) 
(w x 0.02), (iv) and (v). For each weight choice we averaged 
over 10 4 different MSTs. 



of the (i) and (ii) MSTs through the following argument: 
Since the MST tends to avoid links with large weights, 
it effectively shuns the hubs for the cases Wij = kikj and 
Wij = k max , utilizing instead, whenever possible, links 
connecting low degree nodes (Fig. |2t). Consequently, 
all the hubs are marginalized and the MST degree dis- 
tribution must have a narrow range. This argument is 
supported by Fig. an d b, where we show examples 
of MSTs for weight choices (i) and (ii) respectively. The 
sizes of the nodes in the figure reflect their degree in the 
original network. It is evident that the majority of the 
hubs are located on the branches (k = 1 degree nodes) of 



the MST (Fig. HJ. This reliance on the small nodes and 
tendency to avoid the hubs forces the MSTs generated 
by method (i) and (ii) to be very similar to each other. 
Indeed, we find that for a given network but weights cre- 
ated by methods (i) and (ii) , 87% of the links in the two 
MSTs are in common. This explains the similar visual 
appearance of the two MSTs (Fig. [2^ and b). 

The second class of MSTs is well represented by weight 
choices (iii)-(vi), resulting in power-law MST degree dis- 
tributions. The similarity between weight schemes (iv)- 
(vi) is emphasized by the fact that their MST degree 
distributions follow a power-law with the same exponent 
7 = 2.4 (Fig. Eb)- Indeed, the links with the low- 
est weights are now connected to the hubs of the original 
network, and the MST grows utilizing these hubs exten- 
sively. Hence, the hubs of the full network experience 
only a slight reduction in their degree and are found at 
the center of the resulting MSTs (Fig. 0Ji), while the 
intermediate-degree nodes sustain large losses of neigh- 
bors and are found at the surface of the network with 
one or two neighbors (Fig. 0J;). 

The distribution of link weights on the MST also dis- 
plays two distinct behaviors, being either power-law (Fig. 

or exponential (Fig. [3JI). It is interesting to note 
that MSTs with exponential degree distribution (Fig. [3^) 
display power-law weight distributions with exponents 
a = 3.1 (case (i)) and a — 3.0 (case (ii)) (Fig. |3t). On 
the other hand, for weight choices (iii)-(v) the degree dis- 
tribution of the MST is power law and the MST weight 
distribution is exponential or stretched exponential (Fig. 
|3Ji). For (vi) — l/k m %n both the degree and the 
weight distribution of the MST are scale free. Finally, if 
the link weights are distributed uniformly and randomly 
the resulting MSTs have a power law degree distribution 
and an exponentially tempered weight distribution |31| . 

In order to investigate the effect of the degree corre- 
lations on the MSTs for weight choices (i)-(vi), we ran- 
domized the weights of the original network by randomly 
selecting pairs of links and exchanging their weights until 
all correlations between weights and the local network 
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FIG. 4: (Color online) Minimum spanning trees of a TV — 10 3 node scale- free network for weight choices (a) kikj, (b) k max 
and (c) k m i n . The size of a node represents its degree in the full network, and the color of a link represents its weight from 
low (black) to high (green). Note that the MST degree distribution is exponential for (a) and (b) and a power law for (c). 
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topology were lost. Invariably, the resulting MSTs were 
scale-free with a degree exponent similar to that of the 
original network, 7 sa 3. The local structure of the MSTs 
are very different, however, with only 52% of the links 
staying the same in a pairwise comparison between MSTs 
with weight choices (i) and (ii), suggesting that the func- 
tional form of the weight distribution is inconsequential 
for the degree distribution of the MSTs. To understand 
this we recall that only the ranking of the link weights and 
not their absolute value matters 0. Therefore, by re- 
moving the correlations between the local network struc- 
ture and weights we effectively map the problem onto 
that of weights being uniformly random. Indeed, the de- 
gree distribution of the MST in this case is also power-law 
with 7 w 3 31]. However, the MST weight distributions 
continue to depend on the weight distribution of the orig- 
inal network. 

Discussion. — As networks play an increasing role in 
the exploration of complex systems, there is an immi- 
nent need to understand the interplay between network 
dynamics and topology. While focusing on the MSTs 
of scale-free networks, our results emphasize the signifi- 
cance of correlations between link-weights and local net- 
work structure. We find that if correlations are present, 
two classes of MSTs exist, following either a power-law 
or an exponential degree distribution. The removal of 
correlations renders the MSTs scale-free, independent of 
the choice of the weight distribution. This result raises 
interesting questions regarding our ability to quantify the 
influence of weights. 

Our findings could serve as a natural starting point 
towards the systematic exploration of weighted networks. 
For example, while we have assumed that the weights are 
static, incorporating their time-dependence may reveal 
novel dynamical rules. Second, we model the weights as 
solely dependent on the topology, potentially overlooking 
correlations among the weights themselves. Uncovering 
the role of such correlations remains a challenge for future 
research. 
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